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ABSTRACT 

We examine the absorption of cosmic microwave background (CMB) photons by formaldehyde (H2CO) over 
cosmic time. The ^-doublet rotational transitions of H2CO become "refrigerated" — their excitation temper- 
atures are driven below the CMB temperature — via collisional pumping by molecular hydrogen (H2). "Anti- 
inverted" H2CO line ratios thus provide an accurate measurement of the H2 density in molecular clouds. Using 
a radiative transfer model, we demonstrate that H2CO centimeter wavelength line excitation and detectability 
are nearly independent of redshift or gas kinetic temperature. Since the H2CO /if-doublet lines absorb CMB 
light, and since the CMB lies behind every galaxy and provides an exceptionally uniform extended illumina- 
tion source, H2CO is a distance-independent, extinction-free molecular gas mass-limited tracer of dense gas in 
galaxies. A Formaldehyde Deep Field could map the history of cosmic star formation in a uniquely unbiased 
fashion and may be possible with large bandwidth wide-field radio interferometers whereby the silhouettes of 
star-forming galaxies would be detected across the epoch of galaxy evolution. We also examine the possibility 
that H2CO lines may provide a standardizable galaxy ruler for cosmology similar to the Sunyaev-ZeFdovich 
effect in galaxy clusters but applicable to much higher redshifts and larger samples. Finally, we explore how 
anti-inverted meterwave H2CO lines in galaxies during the peak of cosmic star formation may contaminate H I 
21 cm tomography of the Epoch of Reionization. 

Subject headings: cosmic background radiation — dark ages, reionization, first stars — galaxies: high-redshift 
— galaxies: ISM — galaxies: star formation — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Centimeter-wave /if-doublet rotational transitions of 
formaldehyde (H2CO; Figure [U can be collisionally "refrig- 
erated" below the cosmic microwave background (CMB) 
temperature such tha t these transitions absorb CMB photons 
(iPalmer et al.l 11969b iTownes & Cheungl 11969b IE vans et all 
119751: IGarrison et all I1975T " The manner in which the 
centimeter lines of H2CO are cooled is equivalent to a maser 
pumping process: a pump drives an over-population of states, 
but in this case they are lower-energy states, and the level 
populations become "anti-inverted" compared to thermal. 
When the line excitation temperature drops below the local 
CMB temperature, CMB photons may be absorbed. Since 
absorption lines are detectable independent of distance, 
and since the CMB illuminates all gas and provides an 
exceptionally uniform illumination that lies behind every 
galaxy, a Formaldehyde Deep Field (FDF) can provide a 
mass-limited census of the dense molecular gas associated 
with star formation. If multiple lines are observed then 
physical molecular gas densities may also be obtained fairly 
inde pendently of the gas k inetic temperature or other factors 
(e.g. jMangum et al.ll2008h . 

In this Letter, we examine the properties of the anti-inverted 
centimeter-wave H2CO ^-doublet lines as a function of phys- 
ical conditions and cosmological redshift in a feasibility study 
for a census of the cosmic gas evolution and star forma- 
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tion history via an FDF. We also discuss the possibility that 
H2CO lines can provide a standardizable ruler for cosmology 
and may contaminate H I 21 cm tomography observations of 
the Epoch of Reionization (EoR). This work examines radia- 
tive transfer mo dels only, but is supported by observations 
of ou r galaxy dMangum & Wootten 1993: iGinsburg et al.l 
1201 ll) . nearby starburst galaxies (iMangum et al.l 12008b . and 
the gravitational lens ga laxy toward B02 18+357 at z = 0.68 
(IZeiger & Darlingll2010l) . 

We assume 
72 km s" 1 Mpc" 



with H 



a no-curvature cosmology 
1 ,fi„, = 0.26, and fi A =0.74. 

2. RADIATIVE TRANSFER 

We use the RADEX one-dimensional non-LTE radia- 
tive transfer code with the escape pro bability deter- 
mined by the large veloc ity gradient (LVG; ISobolevI 119601: 
iGoldreich & Kwanlll974l) approximation to examine the ex- 
citati on and thus detectability and interpretation of cm H2CO 
lines ( I Van der Tak et al.l 120071). We employ the H2 c ollision 
rate coefficients calculated bv lTroscompt etal.l(l2009bl )forthe 
10 lowest levels of ortho-H2CO (to 7 = 5; FigureHJ for gas ki- 
netic temperatures of 5-100 K. We have compar ed the results 
from this code with another LVG treatment dHenkel et al.l 
119801) an d with RADEX (using the LVG approximation), both 
using the lGreenl (119911) collision rates scaled from He to H2, 
and all are in good agreement for low H2 ortho-to-para ratio 
(OPR). The latter models include the first 40 energy levels of 
ortho-H2CO; agreement between these models a nd those us- 
ing only the 10 levels of iTroscompt et all (l2009bl) lends cred- 
ibility to our treatment of warm molecular gas in star-forming 
galaxies. 

For the radiative transfer calculations we assume statistical 
equilibrium but allow all line excitation temperatures to float 
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FIG. 1. — Formaldehyde energy levels for the first three rotation ladders. 
Arrows indicate allowed transitions, and line frequencies are listed in GHz. 
The /T-doublet levels splitting J states of the K a = 1 ortho-Lb CO rotation 
ladder are the focus of this Letter. The wavelengths of these transitions are, 
from lowest energy to highest, roughly 6, 2, 1, 0.6, and 0.4 cm. 

because the dominant collision partner, H2, "pumps" H2CO 
into non-thermal excitation state distributions. We assume 
that the CMB is the dominant continuum source and that there 
is no significant contribution to the radio continuum from the 
host galaxy itself. An examination of this effect in the star- 
burst galaxy M81 shows a small effect that d oes not quali- 
tatively alter our results dMangum et al.ll2008l) . We also as- 
sume an H2CO/H2 abundance of 10~ 9 , which is reasonable in 
molecular gas-rich galaxies, even at high redshift (Zeiger & 
Darling 2010 obtain an abundance 0. 1-2.5 x 10~ 9 in the grav- 
itational lens galaxy toward B02 18+357 at z = 0.68), and an 

H 2 COOPRof3. 

Troscompt et al.l (l2009allbl) show that the anti-inversion of 



ortho-H2CO is principally driven by collisions with para-H2. 
Since the H2CO line excitation temperatures critically deter- 
mine the detectability and interpretation of lines, the effect 
of the unknown H2 OPR is a concern. To address this, we 
examine model runs with H2 OPR of 3 and 10~ 2 (effectively 
zero) and find small differences in line excitation temperature 
(< 0.3 K) and no significant effect on the broad results and 
interpretation (Section |3}- The basic phenomenon of anti- 
inversion of cm H2CO lines and their ability to produce sil- 
houettes of star-forming galaxies against the CMB is robust 
to the unknown H2 OPR, as is the nearly temperature- and 
redshift-independent nature of the line detectability and exci- 
tation. 

We examine the model predictions for cm H2CO line exci- 
tation for varying n(H2), gas kinetic temperature, and redshift. 
When not varying, we assume w(H2) = 10 4 cm" 3 , Tx m = 40 K, 



and H ? OP R 10" 2 (e.g..lMaret & Berginl (120071) , iPagani et al.l 
( 120091) . and lTroscompt et al l (l2009aF infer ratios in the range 
0.01-0.1 in cold Galactic pre-stellar molecular gas). We fo- 
cus on the first five /if-doublet lines because these produce 
the strongest signals (millimeter lines thermalize in typical 
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FIG. 2. — Formaldehyde line excitation vs. Fb number density. We plot sev- 
eral redshifts to demonstrate the excitation redshift dependence. Line style 
indicates transition, and color indicates redshift. Bottom: excitation temper- 
ature vs. log nfFL;). Middle: excitation temperature normalized by the 6 cm 
line excitation temperature vs. nfFL). We plot all five lines at z. = and the 
2 cm line at redshifts z = 0—4, normalized to the 6 cm line at each redshift. 
Top: observer-frame temperature decrement, A7bb s (EquationfT) for an op- 
tically thick line vs. nCFbJ, We plot redshifts z = 0^1 for the representative 
6 cm and the 1 cm lines (other lines are omitted for clarity). A7ob s = K 
represents the transition from absorption to emission. 



conditions), and the iTroscompt et all (1200 9b) calculations do 
not include higher states. These lines have wavelengths of 
roughly 6, 2, 1, 0.6, and 0.4 cm (4.8, 14.5, 29.0, 48.3, and 
72.4 GHz, respectively). 

The observed line temperature depends on the decrement 
between the rest-frame excitation temperature and CMB tem- 
perature, redshifted to z = 0, and the line optical depth: 



. T x (z) - Tcmb(z) _ t 

A7b bs = : (1-e ) 

1 +z 



(1) 



where T C mb(z) = 7cmb,o(1 +z) and r CM B,o = 2.73 K. For a non- 
unity covering fa ctor, the optical depth would be correspond- 
ingly diminished. iMangum et al J ( 120081) estimate the apparent 
line optical depths in nearby star-forming galaxies: 0.002- 
0.007 in the 6 cm line and 0.0015-0.005 in the 2 cm line. 
IZeiger & Darlingl (120101) estimate 0.017 at 6 cm and 0.008 at 
2 cm in the molecular absorber toward B02 18+357 at z = 0.67. 
The Galactic range among ultracompact H II regions is sim- 
ilar, but with a hi gher upper bound: 0. 006-1.0 at 6 cm and 
0.003-0.3 at 2 cm dGinsburg et al.ll201 ll) . 

The goal of this work is to examine the observer-frame tem- 
perature decrement ( T x (z) - 7cmbCz)) /( 1 +z) of cm H2CO lines, 
which is equivalent to A7bb s for optically thick lines. We fo- 
cus on the excitation temperature of the cm H2CO lines as a 
function of gas kinetic temperature, H2 density, and redshift 
because these factors will most significantly determine the 
utility and interpretation of H2CO line observations at high 
redshift. 

3. RESULTS 
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FIG. 3. — Formaldehyde line excitation vs. redshift. Line style indicates 
transition, and color indicates H2 OPR. Bottom: excitation temperature vs. 
redshift. The CMB temperature is bold. T x < Tcmb indicates absorption of 
CMBphotons. Top: observer-frame temperature decrement ATobs (Equa- 
tion [T) for an optically thick line vs. redshift. The fairly flat trend in these 
decrements shows the nearly redshift-independent detectability of cm H2CO 
lines. 

Our radiative transfer models show that the excitation tem- 
perature of cm H2CO lines is a strong function of H2 number 
density, collisionally driven below 7cmb over roughly three 
orders of magnitude, 10 2 cm" 3 < n(H2) < 10 s cm" 3 (Figure|2]i. 
At low densities the H2CO transitions thermalize to the CMB 
temperature, T x — > Tcmb, because the molecule becomes de- 
coupled from the gas and assumes the radiation temperature. 
At high densities the H2CO excitation thermalizes to the ki- 
netic temperature, T x — > T^ n . 

The range over which n(H2) causes anti-inversion of the 
H2CO lines shifts slightly to higher densities with increasing 
redshift, by ~0.5 dex from z = to z = 4, but this is entirely 
due to the 7cmb = (1 + z) 2.73 K scaling rather than a change 
in physical pumping conditions. The minimum A7bb s like- 
wise grows slightly (becomes more negative) with increasing 
redshift, decreasing by ~0.2 K from z = to z > 1 . 

For fixed n(H2) and T^ n , T x grows roughly linearly with 
redshift (Figure|3). But since the CMB temperature scales as 
Tcmb = (l+z)2.73 K, the redshift dependence largely drops 
out of A7bbs, and the observable line decrement changes very 
weakly with redshift. Thus, not only is there no distance dim- 
ming of anti-inverted cm H2CO lines thanks to the absorp- 
tion of CMB photons, but the excitation conditions are not 
changed significantly by the growth of CMB temperature with 
redshift. 

At each redshift, T^ n has a T^ n = Tcmb = T x floor: the anti- 
inversion disappears when the radiation, excitation, and gas 
kinetic temperatures are equal (Figure 2}. The lower panel of 
Figure |4] indicates a nearly redshift-proportional line excita- 
tion: T x becomes insensitive to T^ n once collisions dominate 
the level populations, and redshift alone dictates T x . Except 
for z ~ 0, where Tqmb is very low, the line decrements are 
insensitive to T^ n after an initial drop from T x = 7kj n = Tcmb- 
In fact, higher energy /if-doublet lines become anti-inverted at 
z > 1. The 6 and 1 cm lines show excellent overlap in A7bb s 
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FIG. 4. — Formaldehyde line excitation versus gas kinetic temperature. We 
plot several redshifts to demonstrate the impact of the (1 +z) CMB temper- 
ature scaling. Line style indicates transition, and color indicates redshift. 
Bottom: excitation temperature versus T K \ n . The black diagonal line shows 
Tx = Tkin, and the convergence point of the transitions at each redshift indi- 
cates where radiation, kinetic, and excitation temperatures are equal. Middle: 
excitation temperature normalized by the 6 cm line excitation temperature 
versus T K \ a . We plot all lines at z, = and the 2 cm line at redshifts z = 0-4 
normalized to the 6 cm line at each redshift. Top: observer-frame temperature 
decrement ATobs (Equation[T} for an optically thick line versus T^ n . We plot 
redshifts z = (Ml for the representative 6 cm and the 1 cm lines (other lines 
omitted for clarity). After the plunge from Tjj,, = T X = Tcmb* T % and ATobs 
depend only weakly on T^ n with the notable exception of the higher energy 
z = lines, which lie on the threshold of anti-inversion. 

for z > 1 and a nearly 7^ -independent trend. 

As the upper panel of Figure |4] shows, there is little change 
in the observed line depth versus redshift or T^ n with the no- 
table exception of the 6 cm line at z — 0, and there is rapid 
convergence of all lines at z > 1 for realistic values of T^\ n 
(> 20 K). Figure H] also demonstrates that there is little vari- 
ation in the excitation temperature ratios between lines for 
z > 1 as a function of gas kinetic temperature. Since line ra- 
tios indicate the H2 molecular gas density, this demonstrates 
that we do not need precise values for kinetic temperature to 
employ the H2CO "densitometer," regardless of redshift and 
regardless of observed lines. 

For nearly all redshifts and in realistic physical conditions, 
the 6 cm line shows the strongest anti-inversion, and the anti- 
inversion monotonically decreases with increasing line en- 
ergy. Remarkably, for most scenarios, the line decrements 
grow from z = to z ~ 1, making the lines more detectable at 
high redshift than locally. 

4. DISCUSSION 

Our model predictions for formaldehyde excitation are 
consistent with observations of molecular clouds in the 
Galaxy, in nearby star-forming galaxies, and in a gravita- 
tional lens. Galactic dark clouds sh ow a 6 cm line d ecre- 
ment of T x - T CM b ^ -0.7 to -0.4 K dEvans et alj|1975l) . and 
an LVG model applied to Galactic giant molecular clouds ob- 
tains comparable li ne decrements, ranging from -1 to K 
dHenkel et al.lll980l) . The H2CO observations and radiative 
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transfer modeling of star-form ing galaxies in the local uni- 
verse by iMangum et al.l ( 120081) also agree with our models. 
Th e molecular cloud studie d toward B02 18+357 at z = 0.68 
by IZeiger & Darlingl (120101) found line decrements of (T x — 
7cmb)/(1 +z) < -2.1 K for the 6 cm line (for 2 x 10 3 cm" 3 < 
«(H 2 ) < 1 x 10 4 cm" 3 ) and -1.5 to -1.8 K for the 2 cm tran- 
sition. 

4.1. A Formaldehyde Deep Field 

An FDF would survey the history of cosmic star formation 
in a fashion similar to other star formation proxies such as 
CO lines, IR/submillimeter luminosity, radio continuum, or 
UV light. But an FDF will differ from traditional proxies 
in three essential ways: (1) line detection is nearly distance- 
independent and thus a survey for H2CO lines would be a 
dense molecular gas mass-limited survey; (2) the H2CO line 
strength is insensitive to gas kinetic temperature and dust tem- 
perature; and (3) H2CO lines are unaffected by dust opacity. 
We suspect that an FDF may be as sensitive to extended gas- 
rich star-forming galaxies such as the "BzK" galaxies as it 
is to compact merger-induced starbursts such as ULIRGs and 
star formation-dominated submm galaxies. While the FDF 
concept has its shortfalls, such as a reliance on radiative trans- 
fer modeling for interpretation, it is a useful survey tool that 
is orthogonal to current star formation survey methods, and it 
overcomes the problem of obtaining redshifts for dusty galax- 
ies identified by their dust continuum. H2CO lines select star- 
forming gas, eliminating the uncertain active galactic nucleus 
dust-heating contribution to continuum star formation prox- 
ies. The strongest signals will arise from the highest molecu- 
lar gas masses, which will include submm galaxies, ULIRGs, 
HyLIRGs, red quasar hosts, dust-obscured galaxies, and BzK 
galaxies. 

The requirements for an FDF are: (1) beam-matched ob- 
servations to the star-forming region(s) in galaxies; (2) large 
bandwidth, which provides a large redshift span; (3) spectral 
line sensitivity; and (4) either observation of more than one 
H2CO line or ancillary data in other parts of the EM spectrum 
to break a one-line redshift degeneracy. Radio interferometry 
meets these criteria and captures many galaxies per primary 
beam, but reduces the brightness temperature sensitivity. One 
should thus work at higher frequencies and use a more com- 
pact array, but higher frequency lines have lower optical depth 
and higher T x (reduced A7bb s ; Figures [2] and [3). A compro- 
mise can be found, likely using the 2 or 1 cm lines, despite 
the intrinsically stronger 6 cm line signal. 

Figure [5] shows the angular size of 2, 5, and 10 kpc molec- 
ular regions and the Jansky Very Large Array (VLA) synthe- 
sized beam for the 1, 2, and 6 cm H2CO lines for redshifts 
0-10. The 4 GHz instantaneous frequency coverage of the C- 
band spans the remarkable redshift range z = 0.8-6.2 — most 
of the epoch of galaxy evolution — because as the 2 cm line 
redshifts out of the C-band, the 1 cm line enters the band- 
pass. Even compact star-forming regions in galaxies could 
be beam-matched in the X-band over a significant redshift 
range for galaxy evolution studies: z = 0.2-0.8 and z = 1.4— 
2.6 simultaneously. A deep observation in one of these bands 
would detect every galaxy above the gas mass limit in the field 
of view and would provide spectroscopic redshifts. 

4.2. Formaldehyde Cosmology 

Formaldehyde silhouettes of galaxies against the CMB of- 
fer a novel probe of cosmology similar to the Sunyaev- 
Zel'dovich (SZ) effect in galaxy clusters. The thermal SZ 
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FIG. 5. — Angular sizes of fiducial 2, 5, and 10 kpc molecular gas regions 
in galaxies (bold black lines) and the VLA A-array angular resolution of the 
1, 2, and 6 cm lines of tbCO versus redshift. Galaxies can be beam-matched 
if the angular resolution track falls below a given physical size locus. Colors 
indicate the VLA wide-band receiver coverage for each line (L-band spans 
1-2 GHz, S spans 2-4 GHz, C spans 4-8 GHz, X spans 8-12 GHz, and Ku 
spans 12-18 GHz). 

effect depends on the electron density as n e while the X-ray 
emission from the cluster gas depends on A comparison of 
the two quantities yields a measure of the line-of-sight length 
scale of the cluster, I. If the proper size and line-of-sight size 
are equal or can be related via models, then the observed an- 
gular size of the cluster yields the angular size (or luminosity) 
distance: 

D A =D L (l+zT 2 =^. (2) 

A similar measurement may be possible using H2CO ab- 
sorption of the CMB in galaxies. H2CO line ratios indicate 
the local H2 density, and line depths yield the H2CO column 
density. Since the column density is simply the space den- 
sity integrated along the line of sight, a comparison of the 
two quantities measures the line of sight physical size of the 
dense gas region of the galaxy. Assuming the line of sight and 
proper sizes to be equal or the galaxy has a measurable incli- 
nation, one can measure Da (D£) from a spatially resolved 
spectral line map. 

Neither clusters nor galaxies are spherical, but many galax- 
ies can be observed per field of view, and in a statistical aggre- 
gate of randomly oriented galaxies, the assumption of equal 
size and length scales will be valid. This redshift-independent 
H2CO technique can employ galaxies far beyond the observed 
Z ~ 1 limit on clusters, providing a unique standard ruler from 
the end of reionization to the formation of galaxy clusters. 

While the number of galaxies with metal-rich gas at high 
redshift may decline compared to z = 2-3, molecules detected 
in quasar hosts up to z = 6.42 have demonstrated that high den- 
sity regions i n the early univers e form metals and molecules 
quickly (e.g.. I Walter et al. 2003). This suggests that one may 
be able to identify standard rulers at very high redshift and 
make unprecedented cosmological measurements over a large 
span of the age of the universe (unlike cosmological redshift, 
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Di and Da are history-integrated quantities). While the star- 
forming regions in galaxies may show size evolution with red- 
shift, this method does not rely on a fixed size for galaxies, 
only on a measurable size. 

The largest possible source of uncertainty for H2CO cos- 
mology is the H2CO abundance, which is needed to convert 
from n(H 2 ) to «(H 2 CO). Since I ~ jV(H 2 CO)/«(H 2 CO), but 
n (H 2 ) = «(H 2 CO)/Xh,co is the quantity derived from line ra- 
tios, the measured size of the molecular region in galaxies is 
inversely proportional to the H 2 CO abundance. 



4.3. Formaldehyde Contamination of H 1 21 cm Tomography 
at High Redshift 

One significant concern for studies of the H I 21 cm line 
before, during, and after the EoR is whether H 2 CO absorp- 
tion of CMB photons may form a contaminating foreground. 
The main lines of interest would be those below or redshifted 
below - 200 MHz (zm > 6): the 2 20 -2 21 71.1 MHz and 
3 2 i -3 22 355.6 MHz lines of para-H 2 CO are the best candi- 
dates (Figure[T|i. The ground state of the K a = 2 rotation ladder 
is 57.6 K above ground, but radiative transiti ons to lower para- 
H?CO states are forbidden. As discussed in lZeiger & Darlingl 
(2010'), warm molecular gas can populate the K a = 2 states 
and collisions with H 2 will drive these meter-wave transi- 
tions to very low excitation temperatures, of order 10-30 mK, 
nearly compensating for the low optical depths of these lines 
compared to the centimeter-wave lines. These low-frequency 
lines have a negligible excitation temperature compared to the 
CMB such that Ar 0bs ^ -(1 -e~ T ) 2.73 K ~ -r • 2.73 K, and 
the observed line depth is simply determined by its optical 
depth regardless of redshift. 

Comparing H 2 CO meter-wave optical depths to cosmolog- 
ical H I 21 cm signals of nomina l brightness temperature 
20 mK (e.g JGnedin & ShaveifcOOl , we find that the 2 20 -2 21 



and 3 2 i -3 22 lines are weaker than the H I 21 cm line for 
t < 0.007. However, our model predicts lines with r ~ 0.2 
and 0.007, respectively. In slightly warmer 60 K gas, the op- 
tical depths rise to 0.4 and 0.02, and both lines become a po- 
tentially problematic source of contamination to H I 21 cm 
tomography. The H I and H 2 CO signals may be separable by 
angular scale and spectral width, but one can imagine unre- 
solved star formation in foreground large-scale structures near 
the peak of cosmic star formation at z — 1-3 creating H 2 CO 
355.6 MHz line signatures that could mimic H 1 21 cm signals 
atz~ 10. 

5. CONCLUSIONS 

Collisions with H 2 cause cm H 2 CO lines to become "anti- 
inverted" to the point where line excitation temperatures 
drop below the rest-frame CMB temperature and thus cre- 
ate silhouettes of the dense molecular gas regions in galax- 
ies against the uniform background. We have demonstrated 
using radiative transfer models that cm H 2 CO lines observed 
in an FDF can provide a distance-independent extinction-free 
mass-limited census of the cosmic history of star formation. 
The detectability and interpretation of H 2 CO line surveys is 
nearly independent of redshift and gas kinetic temperature. 
We have also examined the possibility that H 2 CO lines can 
be used to standardize galaxies as cosmic rulers and that anti- 
inverted meterwave H 2 CO lines may confuse or contaminate 
H I 21 cm tomography of the EoR. Significant new observa- 
tions of H 2 CO at high redshift will soon be possible using 
wide-band interferometric facilities such as the VLA. 

We are indebted to K. Eggert for writing support. We ac- 
knowledge the support of the NSF through award GSSP07- 
0015 from the NRAO and grant AST-0707713. Th e LUNAR 
consortium Qhttp : / / lunar . Colorado . edu) l is funded 
by the NASA Lunar Science Institute (Cooperative Agree- 
ment NNA09DB30A). 
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